Neutrino emission from high-energy component gamma-ray bursts 
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ABSTRACT 

Gamma-ray bursts have the potential to produce the particle energies (up to 10^^ eV) and 
the energy budget (10^'* ergyr"-'^ Mpc~'^) to accommodate the spectrum of the highest energy 
cosmic rays; on the other hand, there is no observational evidence yet that they accelerate 
hadrons. Fermi recently observed two bursts that exhibit a power-law high-energy extension of 
the typical (Band) spectrum that extends to ~ 30 GeV. On the basis of fireball phenomenology 
we argue that they, along with GRB941017 observed by EGRET in 1994, show indirect evidence 
for considerable baryon loading. Since the detection of neutrinos is the only unambiguous way 
to establish that GRBs accelerate cosmic rays, we use two methods to estimate the neutrino 
flux produced when the baryons interact with fireball photons to produce charged pions and 
neutrinos. While the number of events expected from the Fermi bursts is small, we conclude 
that an event like GRB941017 will be detected by IceCube if gamma-ray bursts are indeed the 
sources of the cosmic rays. 

Subject headings: gamma rays: bursts — gamma rays: observations — neutrinos 



1. Introduction 

The sources of the extragalactic cosmic rays 
with energies in excess of '--^SxlO^^eV remain a 
mystery, but one of the best motivated candi- 
dates is gamma-ray bursts (GRBs). Large cosmic- 
ray energies can be achieved in the prompt phase 
of the GRB fireball where internal shocks have 
the potential to accelerate charged particles up 
to ^lO^ieV (iVietri 1995, Waxman 199^. Ad- 
ditionally, the total energy density in the Uni- 
verse of cosmic rays must be matched by a suf- 
ficiently high hadronic energy density in the GRB 
with Phad, GRB ~ Pen- GRB observations iden- 



tify synchrotron photons produced by the elec- 
trons accelerated in the fireball with an energy 
CeETOT ~ 10^^ ergs, where Etot is the total energy 
released by the burst. GRB fireballs also carry en- 
ergy cbEtot in the form of magnetic fields, and, if 
they are the sources of cosmic rays, energy epExoT 
in protons. Assuming equipartition, £b — Se and 
Ee + Cp + es 1- 

GRBs emerge as credible sources for the ul- 
tra high-energy cosmic rays because their observed 
flux can be accommodated with an energy density 
in protons that is similar to that in electrons, or 
£„ Eg. Recent estimates of the local rate of 



GRBs yield a maximum of fig^l Gpc~^ yr"'^ as- 
suming that GRBs follow the star formation rate. 
For a stronger evolution with redshift, the local 
rate can be as low as 0.05 Gpc^'^ yi'^^- Using 
this result we estimate the electromagnetic energy 
density from GRBs to be in the range Pom.GRB ~ 
noCeETOT = 5 X lO'^^-lO'''' ergsMpc^^yr"!. In 
order for GRBs to be the sources of cosmic rays, 
their hadronic energy density Phad,GRB needs to 
produce the observed cosmic-ray energy density 
Phad.GRB = hoepE'TOT ~ 10''^ergsMpc~^yr-^ 
We therefore conclude that Cp/ee « 1 - 20. The 
extragalactic spectrum is actually expected to ex- 
tend to energies below the knee, but remains un- 
observed due to the l arger co i itribu tion of galactic 
cosmic rays, see e.g. lAhlers (l2005l) . If this is the 
case, the average fraction of proton to electron en- 
ergy needs to be larger. 

Fermi recently observed two bursts, GRB090510 
and GRB090902b, that show a statistically signifi- 
cant deviation from the typic al GRB spectrum de- 
scrib ed by the Band function ( Abdo 20091 : Bissaldi 
20091 ) . A flux of high energy events is detected that 
extends to energies '--^30 GeV following a power- 
law spectrum (Table [Ij . This is similar to the po- 
tentially much more luminous burst observed by 
EGRET in 1994, GRB941017 (Gonzalez 2003). In 
this paper we will investigate the possibility that 
the high energy component results from Tr'^-decays 
and is therefore a signature of proton acceleration. 

The main contribution to the initial opacity of 
the fireball comes from the annihilation of photons 
into pairs. The Fermi observation of a non- 
thermal spectrum up to an energy -Emax of tens of 
GeV can be used to constrain the minimum bulk 
Lorentz factor Fm™ required to make the source 
optically thin at the time of the gamma-ray dis- 
play. For all photons with energy E < E^ax the 
condition Ty-y(E) < 1 must be fulfilled where r-y^ 
is the opacity. The observation of photons with 
energies of tens of GeV requires highly relativis- 
tic outflows with F ~ 10'^. Because, on the other 
hand, the observed energy fiux of order 10~*erg/s 
is typical of an average burst, the large boost fac- 
tor implies that the photon density in the rest 
frame of the burst is low. This is a strong ef- 
fect as the photon density is suppressed by F^"'. 
From T^^{E) = 1, we can determine eg by finding 
the electromagnetic energy over the volume of the 



fireball as a fraction of the total GRB energy: 
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where F is the bulk Lorentz factor of the fireball, 
6t is the variability timescale, Ej is the charac- 
teristic gamma-ray energy (which we take to be 
the peak energy of the event), Tgo is the duration 
of the burst, and Etot is the total energy of the 
burst. We therefore estimate for GRB090510 
to be ~ 0.05 and for GRB090902b to be ~ 0.02. 
From the low values of Ce thus obtained we must 
conclude that protons dominate the fireball in or- 
der to accommodate the total energy of a solar 
mass generated by the primary engine. 

In this letter we will first discuss the properties 
of the bursts. We subsequently compute the neu- 
trino flux inevitably produced when the protons 
interact with fireball photons. Their observation 
would provide incontrovertible evidence for the pi- 
onic origin of the additional high energy compo- 
nent in the burst and support the speculation that 
GRB are the sources of the highest energy cosmic 



Table 1 : Spectral parameters of Fermi GRBs with 
power-law components. The total fluence in 
gamma rays F^°'^ is the sum of the Band flu- 



The 



ence and the power-law fluence F^^ 
fluence for GRB900510 is calculated over the en- 
ergy range 10 keV-30 GeV, and the fluence for 
GRB09092b is calculated over the range 10 keV- 
10 GeV. The parameters for GRB941017 are tab- 
ulated in lGonzalea (|2003r ). 
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rays. 

Is a kilometer-scale neutrino telescope such as 
IceCube sensitive enough to shed light on these 
questions? High energy neutrinos are produced in 
the fireball when protons produce pions in inter- 
actions with the low-energy photon field: 



PI 



■ A+ 



n vr"*" 1/3 of the cases 
p-K^ 2/3 of the cases 



(2) 



27, and the 
e+ i^e v^. 



The neutral pions decay as tt*^ 
charged pions decay as 7r+ ^+ 
Here, a single neutrino carries approximately 
l/4th of the 7r+ energy and a photon carries 1/2 
of th e 7r° energy. Th e calculation of the neutrino 
flux ( Waxman 1998h has b een performed in de- 
tail f or the EGRET bursts (jGuettal l20oi iBecker 
20061 ) with the following results: whereas an aver- 



age burst produces only ^ 10~^ neutrinos, bursts 
that are unusually energetic or nearby produce 
an observable flux in a kilometer-scale neutrino 
telescope of order 10 events per year. We suggest 
that the power-law high-energy spectral feature 
identifies such bursts. 

We will compute the neutrino fluxes expected 
in IceCube using two methods: the standard fire- 
ball model, and the bolometric method which re- 
lates the energy in neutrinos from the decay of 
charged pions to the observed photon energy as- 
suming that it is of pionic origin. We will con- 
clude that the neutrino rates from all 3 bursts 
are predicted to be larger than the average, but 
that a burst like GRB941017 extending to GeV 
energy will be observed by IceCube. We remind 
the reader that IceCube observes cosmic neutri- 
nos in a background of neutrinos produced in the 
atmosphere. Given that neutrinos of GRB origin 
are relatively energetic and that the direction and 
time of the events can be correlated to satellite 
alerts, the atmospheric background is suppressed 
and a single neutrino may represent a conclusive 
observation. 

2. Detection of GRBs at high energies 

In fireball phenomenology the main contribu- 
tion to the opacity comes from the annihilation 
of pairs of photons into pairs. The observa- 
tion of a non-thermal spectrum up to -Emax can 
be used to constrain the minimum bulk Lorentz 
factor Fmin required to make the source optically 



thin. From the condition t^7(£'max) < 1 a lower 
limit for T can be derived that requires a mini- 
mum value of the jet boost factor proportional to 
the highest energy £^max observed: 



A{z,5t,(3) 



n(/3). 



(3) 



Here (3 is the spectral index of the high energy 
part of the Band function and St the variability 
time of the GRB. Here w e use the e xpres sions for 
A{z,5t,p) and n{l3) from iLithwic"^ (|200l[ ). 

The first clear detection of a GRB with a power- 
law component in addition to the standard Band 
spectral form was made by the EGRET satellite 
with the observation of GRB941017 in photons 
with energies up to 200 MeV. The observation be- 
ing statistics limited, it is of particular importance 
that EGRET did not detect a cutoff of the power- 
law component and the flu x therefore p otentially 
extended to higher energy ( Alvarez-Mu fiiz 2004|). 
We plot Fmin vcrsus £^MAx for the measured time 
intervals numbered 2, 3, 4 (Figure [T]), described 
by the parameters of Eq. [3l Here A — 10.2166, 
n = 0.211 for bin 2; A ^ 6.6721, n = 0.226 
for bin 3; ^ = 3.3587, n = 0.258 for bin 4, and 
A = 6.7996, n = 0.231 for bin 5. A typical boost 
factor of F ~ 300 is compatible with an extension 
of the high-energy power law to a maximum of 
18.9 TeV. 




[GeV] 



Fig. 1. — -Emax Vs. Fmin for time bins 2-5 of 
GRB941017. Emax = 18.9 TeV corresponds to a 
boost factor of ~ 300, as indicated by the dashed 
lines. Such an energetic photon would exceed the 
Fermi LAT detection capability (--300 GeV). 
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The identification of GRB9410 17 as a poten- 
tially powerful neutrino emitter (|Alvarez-Muiiiz 
20041 ) suggests a comparison with the similar 
bursts observed by Fermi. Time integrated spec- 
tra of the recently observed GRB090510 and 
GRB090902b show statistically significant devi- 
ations from the Band function fit that extend to 
^30 GeV. In Figure [21 a comparison between the 
fluences of the Fermi bursts and the fiuence of 
GRB941017 is presented revealing that the Fermi 
bursts are significantly less bright, especially if 
GRB941017 reached higher values of i?MAx- 
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Fig. 2. — Fiuence of the high energy component 
above 30KeV as a function of the maximum en- 
ergy i?MAx for GRB941017 and observed fluences 
for GRB090510 and GRBG9G9G2b. 



tional fireball phenomenology. Th e calculation is 
rout ine and i s desc ribed in detail in lGuetta ( 2004h 
and Abbasi ( 2009^ . Using th e obs e rved spectra l 
par ameters li sted in lGonzdlezI (l2003l) . lAbdof(|2009l) 
and iBissaldil (|2009l ) the fireball model (Table [U 
predicts 9.7 x 10"^ neutrinos from GRB941017, 
3.5 X 10"'' neutrinos from GRB09G51G, and 1.3 x 
10-2 neutrinos from GRB090902b in IceCube, 
us ing the IC8G neutrino effective area presented 
in Gonzalez-Garcia 



Table 2: Parameters for fir e ball n eutrino fiuxes. 

Abbasil (|2GG9[ ) is defined as 
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Deviations from the Band-only fit in the spec- 
tra are particularly interesting in the context of 
hadronic acceleration within the fireball and re- 
late closely to predictions of neutrino fluxes de- 
tectable on Earth with km^-scale telescopes. If the 
extra high-energy components originate from pi- 
onic photons they provide an optimal benchmark 
for testing models of hadronic acceleration in GRB 
engines. 

3. Neutrinos 

In the hadronic flreball a burst of high-energy 
neutrinos is expected to accompany the prompt 
flux of gamma rays. Assuming that electrons and 
protons are shock-accelerated in the same region, 
the neutrino spectrum can be calculated from the 
observed photon Band spectrum using conven- 



We also expect that if the observed power-law 
components of the photon flux result from the de- 
cay of neutral pions photoproduced by protons 
on fireball photons, then there should also be a 
related flux of neutrinos due to the production 
and decay of charged pions. However, the ob- 
served power-law spectral components are quite 
flat (E~^ — E~^-^) and do not follow the approxi- 
mate E^'^ spectral dependence expected for shock 
acceleration. This spectral behavior is compatible 
with synchrotron radiation from highly relativis- 
tic electrons and positrons, which are produced 
when high-energy photons scatter in the fireball 
photon field. While the spectral behavior is mod- 
ified in this cascading process, the total energy 
contained in the observed photons corresponds to 
the total energy originally going into high-energy 
photon production. Hence, we can assume that 
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the fluence of gamma rays (= f E (j){E) dE ) 
remains constant. Moreover, since both neutri- 
nos and gamma rays originate in pions, we know 
that their energy fluences must be proportional 
oc Fj. If we assume that each gamma ray 
takes one half of the pion energy and each neu- 
trino takes one quarter, and that since there are 
two neutral pions per charged pion there are four 
gamma rays per neutrino (after oscillations), we 
find the pro portionality constant to be approxi- 
mately 1/8 (| Alvarez- Mufii3 12OO2I ). Therefore, if 
we recalculate the gamma-ray fluences from the 
power-law spectral components for different val- 
ues of Emax, we can equate those fluences to the 
neutrino fluence and hence, given the spectrum of 
neutrinos from pion decay, calculate the normal- 
ization of the neutrino flux. Using the IceCube ef- 
fective area, we find the total number of neutrinos 
detected in IC80 from each burst as a function of 
-E'max (Figure[3]). We call attention to the fact that 
a burst with the same parameters as GRB941017, 
extending to E'max > 2 GeV, will produce > 1 neu- 
trino event in IceCube. Given the absence of back- 
ground events over such a short time interval, a 
single event would be statistically significant and 
would constitute discovery of proton acceleration 
in GRBs. 
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Fig. 3. — Neutrino events associated with the 
power-law spectral component from each burst as 
a function of E^ax- 



4. Summary and Conclusions 

Three gamma-ray bursts with statistically sig- 
nificant high-energy power-law spectral compo- 



nents have been detected thus far. The observed 
characteristics of the bursts indicate (Eq. [1]) that 
it is likely that a large fraction of their total energy 
is hadronic rather than electromagnetic. As a re- 
sult, we assume that these bursts are accelerating 
protons and calculate the associated fluxes of neu- 
trinos. While the event rates of neutrinos derived 
from standard fireball phenomenology are small, 
we take the existence of the high-energy power- 
law components as indicative of the decay of 7r°- 
mesons. This allows us to calculate the magnitude 
of the neutrino flux from the related charged pions, 
and we find that a burst like GRB941017 will be 
observable in IceCube if its power-law component 
extends to energies in excess of ~ 1 GeV. 
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